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Abstract 

We study the crossover between liquid and solid electron phases in a two-dimensional har- 
monic trap as the density is progressively diluted. We infer the formation of geometrically 
ordered phases from charge distributions and pair correlation functions obtained via a large 
scale configuration interaction calculation. 
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1 Introduction 



Semiconductor quantum dots [1,2,3] (QDs) are nanometer sized region of space 
where free carriers are confined by electrostatic fields. Almost all QD-based appli- 
cations rely on electronic correlation effects which are prominent in these systems. 
The dominance of interaction in QDs is evident from the multitude of strongly cor- 
related few-electron states measured or predicted under different regimes: Fermi 
liquid, "Wigner molecule" (the precursor of Wigner crystal in two-dimensional 
bulk), charge and spin density wave, incompressible state reminescent of fractional 
quantum Hall effect in two dimensions [2] . 

The few-body problem has been faced mainly via Hartree-Fock method, Den- 
sity Functional Theory, Configuration Interaction (CI), and Quantum Monte Carlo 
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(QMC) [1,2]. Contrary to mean-field methods, CI and QMC methods allow for the 
treatment of Coulomb correlation with arbitrary numerical precision, and, there- 
fore, represent the natural choice for strongly correlated regimes. In addition, CI is 
a straightforward approach with respect to QMC, and gives access to both ground 
and excited states at the same time, its main limitation being its computational cost. 
Here we apply the CI method to the the strongly correlated regime par excellence, 
namely the crossover region between Fermi liquid and Wigner crystal [2]. We find 
evidence of the formation of Wigner molecules, and confirm previous QMC calcu- 
lations [4,5]. 



2 The crossover between liquid and solid phases 



We consider a few electrons in a two-dimensional harmonic trap, which was proven 
to be an excellent model for different experimental setups [2]. The QD effective- 
mass Hamiltonian is: 
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Here m* and k are the effective mass and the static relative dielectric constant of 
the host semiconductor, respectively, Pi = (xj, yi) is the position of the i-th elec- 
tron, lu is the characteristic frequency of the trap. The cylindrical spatial symmetry 
group of the system is D^, and the single-particle QD wavefunctions, the so called 
Fock-Darwin (FD) orbitals [1], have an azimuthal and a radial quantum number, re- 
spectively. The good global quantum numbers of the few-body system are the total 
angular momentum in the direction perpendicular to the plane, M, the total spin, 
S, and its projection along the z-axis, S z . 

We solve numerically the few-body problem of Eq. (1), for the ground state at 
different numbers of electrons, N, by means of the CI method [6]: we expand the 
many-body wavefunction in a series of Slater determinants built by filling in the 
FD orbitals with N electrons, and consistently with global symmetry constraints. 
On this basis, the Hamiltonian (1) is a large, sparse matrix that we diagonalize by 
means of a Lanczos parallel routine [7]. Note that, before diagonalization, we are 
able to build separate sectors of the Fock space corresponding to different values of 
(M, S,S Z ). To this aim, we rely on tabulated Clebsh-Gordon coefficients obtained, 
once for all, via numerical diagonalization of the S 2 matrix (cf. Sec. 2.3 of Ref. [8]). 
We use, as a single-particle basis, up to 36 FD orbitals, and we need to diagonalize 
matrices of linear dimensions up to ks 10 6 . 

The way of tuning the strength of Coulomb correlation in QDs is to dilute elec- 
tron density. At low enough densities, electrons pass from a "liquid" phase, where 
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Fig. 1. Normalized total density of the five-electron ground state as a function of the ra- 
dius, p, for different values of the correlation strength parameter, A. The length unit is the 
characteristic dot radius, ^qd- 

low-energy motion is equally controlled by kinetic and Coulomb energy, to a "crys- 
tallized" phase, reminescent of the Wigner crystal in the bulk, where electrons are 
localized in space and arrange themselves in a geometrically ordered configuration 
such that electrostatic repulsion is minimized [2,3]. The typical QD lateral exten- 
sion is given by the characteristic dot radius £qb = (h/m*uJo) 1 / 2 , £qd being the 
mean square root of p on the FD lowest-energy level. As we keep N fixed and in- 
crease £q D , the Coulomb-to-kinetic energy ratio A = ^qd/«b [°b = (?™*e 2 ) is 
the effective Bohr radius of the dot] [4] increases as well, driving the system into 
the Wigner regime. 

Figure 1 shows the total electron density vs. p, for N = 5 and for different values 
of A. The density is normalized to N, and the ground state remains (M, S, S z ) = 
(1, 1/2, 1/2) across the liquid-solid transition. While in the non-interacting case 
(A — > 0) the charge is homogeneously spread across the dot, approximately be- 
tween p = and p = 2£q D , for A > 4 an outer ring develops while the weight 
in the inner region is stronlgy depleted. For A > 8 a well resolved peak forms at 
p ps 2.5£qd, corresponding to the "freezing" of the electron liquid into a Wigner 
molecule made of five electrons sitting at the corners of a regular pentagon. This 
interpretation is supported by our analysis of the pair correlation function [3] and 
by the excellent agreement between the computed values of ground-state energies 
of our calculation and those obtained via QMC [4,5]. 

The "crystallization" of the five-electron molecule is confirmed by the analysis of 
the spin-resolved density, normalized to the number of electrons with a given spin 
[Fig. 2]. In the liquid phase [left panel of Fig. 2, A = 2], since the number of 
electron is odd, densities for the two spin species of course differ. In the Wigner 
limit (right panel, A = 10), instead, the system becomes classical and the spin 
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radius p (Z QD ) radius p (/ QD ) 

Fig. 2. Normalized spin-resolved density of the five-electron ground state as a function of 
the radius, p. Left panel: A = 2. Right panel: A = 10. The length unit is the characteristic 
dot radius, £qd- 

degree of freedom irrelevant: densites for electrons with spin up or down coincide. 
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